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Abstract 

Using either quasi-static approximation or exact Mie expansion, we characterize the localized 
surface plasmons supported by a metallic spherical nanoparticle. We estimate the quality factor 
Q n and define the effective volume V n of the n th mode in a such a way that coupling strength 
with a neighbouring dipolar emitter is proportional to the ratio Q n /Vn (Purcell factor). The 
role of Joule losses, far-field scattering and mode confinement in the coupling mechanism are 
introduced and discussed with simple physical understanding, with particular attention paid to 
energy conservation. 
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I. INTRODUCTION 

Metallic nanoparticles support localized surface plasmon-polaritons (SPP) strongly con- 
fined at the metal surface ensuring efficient electromagnetic coupling with neighbouring 
materials, offering a variety of applications such as surface-enhanced spectroscopies^, 
photochemistry^ or optical nanoantennas^. This also opens the way towards control of 
light emission at the nanoscale£~— . The ratio Q/V e ff is generally used to quantify the cou- 
pling strength between a dipolar emitter and a cavity mode (polariton). Q and V e ff refers 
to the mode quality factor and effective volume, respectively. High ratio Q/V e ff may lead 
to strong coupling regime with characteristics Rabi oscillations revealing cycles of energy 
exchange between the emitter and the cavity In the weak coupling regime, the emitter 
energy dissipation into the cavity mode is non reversible and follows the Purcell factor: 



where 7 is the emitter spontaneous decay rate into the cavity compared to its free-space value 
7o, n B is the optical index inside the cavity and A is the emission wavelength. Spontaneous 
emission rate in complex systems follows the more general Fermi's golden rule that expresses 
7 as a function of the density of electromagnetic mode^. However, for describing the 
coupling to a cavity mode, the Purcell factor is usually preferred since it clearly introduces 
the cavity resonance quality factor Q and the mode extension V e ff on which coupling remains 
efficient, bringing therefore a clear physical understanding of the coupling process. 

In this context, we propose to determine the quality factors and effective volumes of 
localized SPPs supported by a metallic nanosphere. Indeed, these quantities are usefull 
parameters to understand and evaluate the coupling mechanisms between dipolar emitters 
and plasmonic nanostructures^ 1 ^. This would help for achieving strong coupling regime' 1 ^ 
or for designing plasmonic nanolasers^ - — . We have chosen a spherical particle, a highly 
symetrical system, since it is fully analytical and simple expressions can be derived with 
clear physical meaning^— . Moreover, our results could be extended to more complex 
structures". 

In section [Til we focus on the dipolar mode and present in details the derivation of its 
quality factor and effective volume. We extend our approach to each mode of the particle 
in section IHH Finally, we discuss the coupling efficiency to one of the particle modes in the 
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last section. For the sake of simplicity, all the analytical expressions are derived assuming a 
particle in air and dipolar emitter perpendicular to the particle surface. The generalisation 
to arbitrary emitter orientation and a background medium of optical index % is provided in 
the appendix. Exact calculations presented in the main text using Mie expansion correctly 
include both the background medium and dipolar orientation. 



II. DIPOLAR MODE 



We first characterize the dipolar mode of a spherical particle. For sphere radius R small 
compared to the excitation wavelength A = 2ttc/u, the electric field is considered uniform 
over the metallic particle. The metallic particle is polarized by the incident electric field Eo 
and behaves CIS db dipole 

pW(w) = 4vre a 1 (a;)Eo, (2) 
e m (uj) — 1 3 

where ol\ is the nanoparticle quasi-static (dipolar) polarisability and e m is the metal dielectric 
constant. The dipole plasmon resonance appears at uj\ such that e m (u>i) + 2 = 0. In 
case of Drude metal, the dipolar resonance is U\ = to p / \f?> with ui p the bulk metal plasma 
angular frequency. However, expression (|3]) does not satisfy the optical theorem (energy 
conservation). It is well-known that this apparent paradox is easily overcome by taking into 
account the finite size of the particle and leads to define the effective polarisability 26 : 

-i 



eff 

a-i 



2k A 



ax ,(k= 2vr/A) . (4) 



The corrective term (2/c 3 ai/3) is the so-called radiative reaction correction and microscop- 
ically originates from the radiation emitted by the charges oscillations induced inside the 
nanoparticle by the excitation field 27 . 

The dipolar polarisability presents a simple shape near the resonance if the metallic 
dielectric constant follows Drude model 20 (co p and T a b s refers to metal plasma frequency and 
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Ohmic loss rate, respectively)^; 

£m = 1 - -^-^ ; (5) 

otf*{u)~- % -R\ (6) 

2(fc 1 /?) 3 u; 1 

Ti = T abs H , (fci = wi/c) . (7) 

where I\ is the decay rate of the particle dipolar mode, and includes both the Joule (T a b s ) 
and radiative [r™ d = 2(/c!.R) 3 a;i/3] losses rates. 



A. Quality factor 

The dipolar response can be described by either the extinction efficiency Q ex t, propor- 
tional to Im(ai), or scattering efficiency Q sca tti proportional to |ai| 2 , and therefore follows a 
Lorentzian profile centered at uj\ and with a full width at half maximum (FWHM) Awi = T\. 

Q«t(w) = ^aft)] (8) 



oc 



«i (w - wo 2 + (ry2) 



2 ' 



Q^M = JLk*\a? f (u)\ 2 (9) 



1 

oc 



(u - Wl ) 2 + (ry2) 2 • 

The quality factor of this resonance is therefore: 

= (10) 
1 i 

As an example, we consider a R = 25 fim silver sphere in air. The Drude model param- 
eters are hu p = 9.1 eV and hT^s = 18 meV that lead to a resonance peak at = 5.2 el/ 
(co>i = 7.98 10 15 Hz) and radiative energy hT r 1 ad = 1 eV. The quality factor is then Qi = 5. 
Note that the metal optical properties are better described when including the bound elec- 
trons in the Drude model: e m = — w 2 /(a; 2 + iT a b s uj), (eoo = 3.7 for silver, see also the 
appendix). In that case, we obtain hu>i = 3.74 eV (ui = 5.7 10 15 Hz), TiT r l ad = 0.13 eV and 
Qi = 24. 

Although Drude model qualitatively explains the shape of the resonance, a more repre- 
sentative value of the quality factor can only be determined using tabulated data for the 




FIG. 1: Extinction efficiency for the dipolar resonance calculated keeping only the dipolar mode 
(n = 1) in the Mie expansion and using experimental value for the metal dielectric constant, a) 
Silver particle in air. b) Gold particle in air. c) Gold particle in PMMA (optical index n# = 1.5). 
The particle radius is R = 25 nm for each case. Fit refers to a Lorentzian fit using parameters 
indicated on the figure. 



dielectric constant of the metal^. The extinction efficiencies associated to the dipolar reso- 
nance are represented in Fig. [TJ For silver particle (Fig. [T^,), the extinction efficiency closely 
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follows a Lorentzian shape, as expected, with a quality factor Qi = 21, in good agreement 
with the value obtained using Drude model (with the contribution of the bound electrons 
included). In case of gold, the resonance profile is not well defined due to interband tran- 
sitions for oo > 4.10 15 Hz but a quality factor of 7 can be estimated (Fig. [Ib). Figure [Tb) 
represents the extinction efficiency for a gold particle embedded in polymethylmethacrylate 
(PMMA) into which metallic particles are routinely dispersed. This leads to a small redshift 
of a resonance, avoiding therefore the resonance disturbance by interband absorption and 
we recover partly the Lorentzian profile 3 -^. 



B. Effective volume 



Coupling rate of a dipolar emitter to a dipolar particle expresses for very short emitter- 
particle distances d (zq — R + d is the distance to the particle center) as^r— : 

for a dipole emitter orientation perpendicular to the nanoparticle surface. Using Eq. EJ 
we obtain, in case of a dipolar emitter emission tuned to the dipolar particle resonance 
(A = Ai = 27tc/Vl) 

it 6^ 3 3 

In order to determine the dipolar mode effective volume, we now identify the coupling rate 
7/7o to the Purcell factor (Eq. [TJ assuming n B = 1), so that we obtain 

vt=*4' (i3) 

in full agreement with expression recently derived by Greffet et al by defining the optical 
impedance of the nanoparticle antenna 31 . For a R = 25 nm radius sphere in air, we estimate 
the mode effective volume = 3.7 10~ 4 /im 3 = (72 nm) 3 for an emitter 10 nm away from 
the particle surface (z = 35 nm). 

Unlike here, usual definition of the mode effective volume does not include the emitter 
position. For instance, in case of cavity quantum electrodynamics (cQED) applications, 
it can be expressed as V e ff = J e|E(r)| 2 c?r/maa;(e|E(r)| 2 ), so that it directly characterizes 
the mode extension. However, in that case, Purcell factor expression (Eq. [T|) is only valid 
for an emitter located at the cavity center where the mode-emitter coupling is maximum 



(mode antinode). In close analogy with the definition used in cQED, we derive later another 
expression for the SPP effective volume, independant on the emitter position (see Eq. [24J, 
and the discussion below). 




III. MULTIPOLAR MODES 

b) 




FIG. 2: a) Dipolar and b) quadrupolar mode profiles of a R = 25 nm silver sphere, embedded in 
PMMA calculated using exact Mie expansion. Silver dielectric constant is taken from ref.— . Color 
and arrows refer to electric field intensity and vector respectively. 



If the excitation field is generated by a dipolar emitter (fluorescent molecules, quantum 
dots, ...), it cannot be considered uniform anymore and the dipolar approximation fails. 
One needs therefore to consider the coupling strength to high order modes (Fig. [2]). Here 
again, we discuss the mode quality factors and volume first using quasi-static approximation 
and then discuss their quantitative behaviour using exact Mie theory. 

The n th multipole tensor moment of the metallic particle is given by: 

p( " , -(2^ a " V "" lE °' < W > 

with (2n + 1)!! = 1 x3x 5 x ... x (2n + 1) and V is the vector differential operator. As 
discussed above, the dipole moment p^ 1 * 1 = 47reoaiE is the unique mode excited in an 



7 



uniform field. However, the dipolar emitter near-field behaves as 1/r 3 and strongly varies 
spatially so that higher modes can be excited in the particle. 

For a Drude metal, the n th resonance appears at oj n = oj p ^Jn/(2n + 1). Therefore, 
higher order modes accumulate near oJoo = oj p /\/2. Moreover, as discussed for the dipolar 
case, quasi-static expression of the n th mode polarisability (Eq. fT5l) does not obey energy 
conservation. Applying the optical theorem, we recently extend the radiative correction to 
all the spherical particles modest. This leads to 



a e J f 



2n+l 1 -1 



(n + l)k 

Oi 

n(2n-l)!!(2n+l)!! 



a n . (16) 



A. Quality factors 

It is now a simple matter to generalize the dipolar mode analysis reported in the previous 
section to all the particle modes. The n th mode polarisablity can be approximated near 
resonance. A simple expression for the polarisability is achieved considering a Drude metal: 

< ff ~ V( ^—-R^ , (17) 

Un 2{0J n — Oj) — 11 n 
n J- abs ~r i n ? 

1 n ~ Un n(2n - l)!!(2n + 1)!! ' { n ~ n/ ] ' 

where T n is the total decay rate of the n th mode, that includes both ohmic losses and radiative 
scattering. As expected, for a given mode n, the radiative scattering rate T™ d oc R 2n+1 
increases with the particle size since it couples more efficiently to the far-field. For instance, 
we obtain hT r ^ ld = 39 meV {UY™ d = 1.8 meV with = 3.7) for the quadrupolar mode of 
a R = 25 nm silver particle in air. As expected, the radiative rate of the quadrupolar is 
strongly reduced compared to the dipolar mode. 

The quality factor associated to the n th mode is therefore 

L n L abs ' L n 

Figure |3] details the quality factor the two first modes of a silver sphere in PMMA using 
tabulated value for e m . The quadrupolar mode presents a quality factor almost 5 times 
higher than the dipolar mode since it has limited radiative losses. Indeed, quadrupolar 
mode poorly couples to the far-field. Finally, similar Q values (Q n ~ 50) are obtained for 
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all the higher modes {n > 3). Here again, assuming a Drude metal and in the quasi-static 
approximation, we qualitatively explains this result. Actually, Q factors of high order modes 



are absorption loss limited (T„ <C Tabs, so called dark modes) and tend to Qoo = ^oo/r 



abs- 



However this overestimates the mode quality factor (Qoo ~ 207) as compared to the value 
deduced using tabulated value for e m and exact Mie expansion. 



a) 



b) 
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FIG. 3: a) Extinction efficiency of a R = 25 nm silver sphere embedded in PMMA, calculated 
using full Mie expansion, b) Contribution of the dipolar (n=I) and quadrupolar (n=2) mode to 
the full extinction efficiency. The quality factor is indicated for each mode. Dielectric constant for 
silver is taken from ref.—. 
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B. Effective volumes 



Last, we express the total coupling strength of a dipolar emitter to the spherical metallic 
particle for very short separation distances (kz <C l)^^: 

Tit ^3 _l_f^ + l) 2 



To 2 (^o) — { z 



(2n+l) 



-Im(agf) . 



(19) 



So that the coupling strength to the n th mode is easily deduced as 

3 1 (n + 1) 2 



1± 

70 ~ 2 {kz f 4 2n+1 ) 



-Im(a 



eff\ 
n J ? 



3 R 



2n+l 



(n + 1) 



<*, 2 (k n Z f JPn+D 



"Qn ■ 



(20) 
(21) 



where we have used approximated expression ffTTj) for the n th polarisability. The mode 
effective volume is then straightforwardly deduced from comparison to the Purcell factor 
(Eq. [TJ) 



v- 

' 1 > 



4nz 2n+4 



(22) 



(n + l) 2 R 2n+l ' 

We obtain for instance quadrupole effective volume V^- = 4ttZq/(9R 5 ) = 3.2 10~ 4 /im 3 = 
(68 nm) 3 for an emitter 10 nm away from a 25 nm radius metallic particle in air. The 
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FIG. 4: Effective volume of the dipole (Vi), quadrupole (V2) and hexapole (V3) modes of a R = 
25 nm sphere as a function of the emitter-particle distance. The emitter is perpendicular to the 
particle surface. 



calculated effective volumes of the three first modes are shown on figure HI At long distance, 
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the smaller effective volume (most efficient coupling to an emitter) is associated to the dipole 
mode since it presents the largest extension. When the emitter-particle distance decreases 
the coupling strength to the quadrupolar (d < 12 nm), then hexapolar (d < 10 nm) mode 
becomes stronger as revealed by their lower effective volume. Most efficient coupling between 
an emitter and one of the particle modes obviously occurs at contact since plasmons mode 
are confined near the particle surface. 

Finally, it is worthwhile to note that mode effective volume is generally defined indepen- 
dently on the particle-emitter distance so that it gives an estimation of the mode extension. 
This is done by evaluating the maximum effective volume available, and for a random emit- 
ter orientation. In case of localized SPP, this is achieved for contact (z Q = R). Since the 
decay rate of a randomly oriented emitter expresses 7 = (7 1 - + 2y")/3, it comes 

11 9 

+ 777777"; > (23) 



V n 3V r HR) 3V»(R) 

V «= (2n + l)(n + l) V °» < M > 

where Vq = 4ttR 3 /3 is the metallic sphere volume and Vn refers to a dipole parallel to the 
sphere surface (see Eq. [36] in appendix, with €b = 1)- Dipolar and quadrupolar mode 
volumes are V\ = 3/2 Vq and V% — 3/5 Vq, respectively. The expression ( l24j) . derived for 
a Drude metal, quantifies an extremely important property of localized SPPs; their effective 
volume does not depend on the wavelength and is of the order of the particle volume. Metallic 
nanoparticles therefore support localized modes of strongly subwavelength extension. High- 
est order modes have negligible extension (V n — > for n — > 00, see Fig. |5]). As a comparison, 
photonics cavity modes are generally limited by the diffraction limit so that their effective 
volume is at best of the order of (X/ub) 3 (see ref.— ). Nevertheless, the extremelly reduced 
SPP volume is achieved at the expense of the mode quality factor. It should be noticed that 
a low quality factor indicates a large cavity resonance FWHM so that emitter-SPP coupling 
can occur on a large spectrum range. 

Note that mode effective volume is generally defined by its energy confinement V e // = 
J e|E(r)| 2 dr/maa;(e|E(r)| 2 )^. Sun and coworkers used this definition and obtained^ V n = 
6Vo/(n + l) 2 that is in agreement with our expression for the dipolar mode volume {V\ = 
1.5Vo) but leads to slightly different values for other modes (e.g V 2 ~ 0.67Vo instead of 
V2 = O.6OV0). Recently, Koenderink showed that defining the mode volume on the basis of 



11 



> 

> 



O 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




1 



2 



3 



4 



5 



N# 



FIG. 5: Normalized effective volume V n /Vo as a function of mode order. 



energy density could lead to understimate the Purcell factor near plasmonic nanostructures^ 4 - 
(however, he defined the coupling rate to the whole system rather than considering the 
coupling into a single mode). Oppositely, we adopt here a phenomenological approach 
where the mode volume is defined so that Purcell factor remains valid. Nevertheless, both 
methods lead to very similar results for the effective volumes of localized SPPs supported 
by a nanosphere. Since mode volume is generally a simple way to qualitatively characterize 
the mode extension, expressions derived here or by Sun et al could be used. 

IV. $ FACTOR 

Purcell factor quantifies the coupling strength between a quantum emitter and a (plas- 
mon) mode but lacks information on the coupling efficiency as compared to all the others 
emitter relaxation channels. For the sake of clarity, it has to be mentionned that coupling 
strength to one SPP mode corresponds to the total emission decay rate induced by this 
coupling. It does not permit to distinguish radiative (jrad) and non radiative (jnr) cou- 
pling into a single mode. This could be done by numerically cancelling the imaginary part 
of the metal dielectric function (see also ref.— for similar discussion in case of coupling to 
delocalized SPPs). However, one can determine the coupling efficiency into a single mode 
as compared to all others modes. This coupling efficiency, or the so-called 0- factor, is easily 
estimated in case of a spherical metallic particle since all the available channels are taken 
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into account in the Mie expansion. Coupling efficiency into n th mode writes 

Pn = , (25) 

< Itot > 

where < 7 >= (7- 1 - + 2^)/3 is the devay rate of a randomly oriented molecule. /3-factor 
is represented on Fig. [6] for the first three modes. A maximum efficiency of 90% can be 
achieved in the dipolar mode (d ~ 10 nm) and 87% into quadrupolar mode (d ~ 15 nm). 
The coupling efficiency into the hexapolar mode is lower (~ 60% around d ~ 15 nm) since 
it has a very low extension, as indicated in Fig. For very short distances, all the coupling 
efficiencies drop down to zero since all the higher order modes accumulate in this region, 
opening numerous alternative decay channels. Moreover, it is possible to efficiently couple 
the emitter to either the dipolar or quadrupolar mode, by matching the emitter and mode 
wavelengths. This is of strong importance when designing a SPASER (or plasmon laser )^ so 
that the active mode can be tuned on the dipolar or quadrupolar mode. This last spasing 
mode would consist of an extremely localized and ultrafast nanosource^. 



CO. 
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FIG. 6: Coupling efficiency into the dipolar (/3i), quadrupolar ((32) and hexapolar (^3) mode 
of a R = 25 nm silver sphere in PMMA, calculated using exact Mie expansion. The emitter 
wavelength is assumed to match the considered mode. Therefore Pi, P2 and /S3 are calculated 
at Ai = 435 nm (dipolar resonance), A2 = 387 nm (quadrupolar resonance), and A3 = 375 nm 
(hexapolar resonance), respectively. 
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V. CONCLUSION 



We explicitly determined the effective volumes for all the SPP modes supported by a 
metallic nanosphere. Their quality factor is also approximated in the quasi-static case or 
calculated using exact Mie expansion. Rather low quality factors ranging from 10 to 100 can 
be achieved, associated to extremely confined effective volume of nanometric dimensions. 
This results in high Q/V e ff ratios indicating efficient coupling strength with a quantum 
emitter. At the opposite to cavity quantum electrodynamics where the coupling strength is 
obtained on diffraction limited volumes thanks to ultra-high Q factors, plasmonics structures 
allow efficient coupling on nanometric scales with reasonable Q factors, defining therefore 
high-bandwidth interaction. Finally, high coupling efficiencies (80% — 90%) to dipolar or 
quadrupolar mode can be achieved and are of great interest for nanolasers realization. 
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VII. APPENDIX 



In this appendix, we derive the expression of mode quality factor and effective volume 
for a spherical particle embedded in homogeneous background of optical index n B = y/e B . 
For a better description of the metal optical properties, we include the contributions of the 
bound electrons into the metal dielectric constant 

^ + UabsU 
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A. Effective polarisability associated to n th SPP mode 

The n th multipole tensor moment of the metallic particle is given by 

rttnV^Eo, (27) 



(n) 47re e B . 



(2n-l)!! 

a n = n{e r~Tl Ri2n+1) > ( 28 ) 

ne m + (n + l)e B 

The resonance angular frequency is then tu n = tu p y / n/[ne 00 + (n + l)es] (woo = 
w p / v^ooH-cb ) • Finally, the effective polarisability, including finite size effects write 22 : 



a eff 



7i(2n-l)!!(2n + l)!! 



a„ • (29) 



with /eg = n^/c the wavenumber in the background medium. 

Considering a Drude metal, we achieve a simple approximated expression for near 
a resonance 

a«" (2n + 1)eg ^ , (30) 

n u>n ntoo + (n + l)e B 2(w n - oj) - iV n 



Tn Tabs ~\~ T n 



2n+l 



r rad (2ri + l)eg , (71 + 1) (fc B fl) 



ne 00 + (n + l)e B n n(2n - l)!!(2n + 1)!! ' 

So that the quality factor expression Q = uj n /T n remains valid but with the corrected 
expressions for resonance frequency u n and total dissipation rate T n . 

B. Effective volumes 

The total coupling strength of a dipolar emitter to the spherical metallic particle, em- 
bedded in riB medium expresses 2 *^: 

\2 



oTK^ES^lfM^), (31) 



Itot 



^Et^M<»). (32) 
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The coupling strength to the n mode is 
In _ 3 1 (n + 1) 



2 



/m(«5) , (33) 



n Bl0 2 (fc B z ) 3 4 2n+1) 



A y i?2 n+i ( 2ri + l)e B (n + 1) 2 



u> n 4vr 2 V^b / 4vr4 n6oo + ( n + l)e B z£ n+1) 
4 ^ 1 ^ + rrW^), (34) 



n Bl0 4(k B z y 4 2 "+D 

/A\ 3 i? 2n+1 (2n + l)e B n(n + 1) 
Bk4 n £oo + (n + l)e B zi> +1) 



<*. 4vr 2 V^b / 8vr4 n6oo + ( n + l)e B z J 2n+1) 
and the mode effective volume is deduced from comparison to the Purcell factor (Eq. [I]) 

= neoo + (n + l)e B Anzl n+i 

(2n + l)e B (n + l) 2 R 2n + 1, K ' 

I, = nepo + {n + l)e B 8tt^" +4 

(2n + l)e B n(n + l)R 2n + l ' 1 ' 

Finally, we define the mode volume as 

1 - 1 I 2 (37 ) 
V n 3^(i?) SKl 1 ^)' 

_ nepo + (n + l)e B 9 TA 

(2n + l)e B (2n + l)(n + l) K °- ^ 

All the expressions obtained in this appendix reduce to the simple analytical case discussed in 
the main text for = 1 (bound electrons contribution neglected) and e B = 1 (background 
medium is air). 
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